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ABSTRACT 

Current flow with density j > nec (c = /r /m. = ion acoustic velocity) 
should encounter resistance due to ion S acoustic turbulence when T >> T.. 

0 i 

This has been observed in experiments with a streaming cesium plasma (B = 0) , 
in which electron current, potential rise due to turbulent resistivity, spec- 
trum of unstable ion acoustic waves, and associated electron heating were 
all measured directly. Kinetic theory calculations for an expanding, un- 
stable plasma, give results in agreement with experiment. The effective 
electron collision frequency is of order 5 to 10 to ^ In a strong magnetic 

field, with T /T. % 1 and current densities less tfian nea (a = /2T / m ), 

01 0000' 
typical for present Tokomaks, the plasma is stable to ion acoustic but unstable 

to current driven electrostatic ion cyclotron waves. Relevant characteristics 

of these waves are calculated and it is shown that for ion 8 > m /m., the 

electromagnetic ion cyclotron wave has a lower instability thresfiolb than the 

electrostatic one. However, when ion acoustic turbulence is present 

(T £ /T. > 5) experiments with Double Plasma (DP) devices show rapid anomalous 

heating of an ion beam streaming through a plasma. These DP devices also 

produce large amplitude, collisionless ion acoustic shocks, both laminar and 

turbulent. At low Mach numbers, a beam of reflected ions is observed. 

Turbulence, either natural or triggered by fluctuations of injected ion waves, 

scatters the reflected ions and heats the shocked ions (observed broadening 

of distribution function). Shock structure is very sensitive to T /T. and 

especially to the presence of light ions: a few per cent He in an A plasma 

can significantly change both the shock amplitude and the laminar/ turbulent 

transition. At higher density ratios (obtained with photoionization and 

microwave, upper hybrid, resonant heating) the Mach number increases (maximum 

observed % 1.8) and the structure washes out. 
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1 . Introduction 

We report here four related experimental and theoretical studies of ion 
acoustic and ion cyclotron turbulence and the role they play in anomalous 
resistivity, viscosity, diffusion and heating and in the structure of colli- 
sionless electrostatic shocks. 


2. Anomalous Resistivity Due to Ion Acoustic Turbulence . 

Experiments have been carried out on a quiescent, streaming, inertially 
confined cesium plasma with T /T . » 1 and no external magnetic field. The 
general configuration and dimensions are shown in Fig. 1 with the plasma 
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Fig. 1. Cesium ion beam configuration (on left) and enlarged view of 
plasma source (on right). 


source [1] shown schematically on the right of Fig. 1. The collector 
electrode, 2 m from the source, may "float", drawing electron and ion cur- 
rents, I Q , of order 6 mA; or it may be biased positively, resulting in an 
electrorj current flow, I, from neutralizer wires to collector electrode, 

I Q < I ^ 10 Iq. For I = Iq, the plasma potential on_axis (measured with 
emissive probes) is nearly constant (dV/dz 'v 10 -10 V/cm) but as I is 

increased, a sharp rise in potential is observed (Fig. 2), at values of z 
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Fig. 2. Anomalous potential rise 
due to electron current flow. 


Fig. 3. Amplitude of ion acoustic 
noise (in arbitrary units) 
measured by Faraday cup. 
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which decrease with increasing I. Faraday cup measurement of ion acoustic 
noise show an abrupt rise in the same region (Fig. 3) and Langmuir probes 
show a correlated increase in electron temperature from the initial 0.25 eV 
to values comparable to AV. Growth of ion acoustic waves launched at the 
accelerating grid is also clearly observable when I > Iq. 

Our theoretical analysis assumes, in conformity with experimental 
measurements, that plasma, density (measured with Faraday cups) and current 
density (measured with small flux coils) are uniform inside a cone, r = 
(Tq/z^Jz, (Zq - 10 cm) so an observer moving with the ion streaming velocity, 
Vq - 2x10 6 cm/sec, sees a density decreasing with time. Because v^ >> c g = 

(T /M) 1 / 2 ^ 6><10 4 cm/sec, there is little propagation from one plasma refion 
to another during the beam transit time and an initial-value formulation 
in the ion rest frame is appropriate, with lab frame frequency coming essen- 
tially from Doppler shift co = k (v^ ± c ) - k^Q. Linear theory predicts a 
wave growth from z^ to z given §y 

z 

log [6E k (z)/E k (z Q ) ] =2 / Im 0 J k dz/V Q = (a/2)(k z /k) h(K,y) (1) 

Z 0 

where <5E, (z ) is the fluctuation electric field at injection; k = k/k (z_) ; 
y = z/z Q ; a = (I/I Q ) (Trm/2M) 1 ' 2 k D (z Q ) - 10-15; k D (z Q ) = (^n(z 0 )e/T e ) 1 ' 2 ; 

h(K,y) = Ky [1 + K 2 y 2 ] l7z - (1 +K 2 ) _l72 (2) 

In the anomalous resistivity region, the experimentally observed growth of 
ion acoustic turbulence (Fig. 2) agrees with Eq. (1) for the wave numbers of 
maximum growth [k - K = y 2 / (1 + y 2 / 3 ) -1 / 2 ] but smaller wave numbers grow 

faster than predicted by Eq. (1), suggesting the presence of strong mode 
coupling as observed in other ion-acoustic wave experiments [2]. The elec- 
tron velocity diffusion coefficient associated with this turbulence, calcu- 
lated in the usual way [3] from the second order Vlasov equation, causes an 
anomalous viscosity for the electrons. In our constant current configuration, 
this must be overcome by a potential rise AV (Fig. 2). If E is the associ- 
ated electric field, the effective collision frequency, defined by V = eE/mV n , 
is given by 

V(z)/oo pi (z 0 ) = (M/m) 1 / 2 G(y)/N D (z Q ) 
where v Q = (I/I Q - l)v Q ; N D (z) = n(z)/K^(z); and 

G(y) = K 5/2 (l + K 2 )" 374 e CtH /16Tr/3a 3/2 H(y 2 -y 2/3 ) 

with H = h(K,y). Typical values of v(z)/w . (z) = v(z)/yw ^( z q) are 5 to 10. 

In deriving (4), we use a saddle point appPoximation valiS^or large y and 
neglect variations in T . We also take < ] 6 E k ( z q ) I 2 > to correspond to bare 
particle discreteness, lince the plasma transit time (y 100 ys) is too short 
for formation of Debye shielding clouds. The early portions of the curves 
in Fig. 2 are in qualitative agreement with (3), but the experimental curves 
rise more rapidly with z, probably due to the aforementioned mode coupling. 

The calculated electron heating is also in agreement with the experimental 
values. In addition, it causes an enhanced beam divergence which appears to 
be the principal cause of the saturation of the anomalous potential rise in 
Fig. 2. 


(3) 

(4) 
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3. 


The Ion Cyclotron Instability 


We report four basic extensions of the Drummond-Rosenbluth [4] investi- 
gation of the electrostatic ion cyclotron instability: arbitrary T /T. = 

electron- to- ion temperature ratio; growth rates well above marginal 6 
stability; finite 3^ corrections, where 3^ is the ratio of ion particle to 
magnetic field pressure; and a comparison of the electrostatic and electro- 
magnetic ion cyclotron current instabilities [20]. 


3.1 Electrostatic Ion Cyclotron Instabilities. 


We have chosen isotropic Maxwellian distributions drifting along a 
uniform-background magnetic field. When T /T. =1, the wave that is marginally 
stable to the smallest current has k 2 R 2 - 1, u) - 1.2A ., and k,/k.. = 12, 
where R. is the thermal ion Larmor radius, w is the wave frequency and k, , 
k )( are wavenumbers perpendicular and parallel to [5], [6], [7]. Here we 
focus on unstable ion cyclotron waves. The dispersion relation [4], [5] 
depends upon k^R? , k A /k'„, V /a [the electron drift relative ions in units 
(2T /m ) 1 ' 2 ] , m/A, T /T . , ana tHe ratio of ion plasma to cyclotron frequency, 
U)2 i^ci* ^ 6 e -*-i m -‘- na ^ e some dependences: we take u) > 10, since this 

implies weak us 2 JQ, 1 ^ dependence [8]. We eliminate Sepenclences upon k^/k,, 
and k^R 2 by selrching numerically for local maxima in the growth rate. 

Choosing M/m = 1836, we plot in Fig. 4 contours of maximum growth rate for 




Fig. 4. Contours of maximum 
growth for the fundamental (n=l) 
ion cyclotron wave (solid) and 
ion acoustic wave (dashed). 


Fig. 5. Maximum growth rate con- 
tours for electromagnetic ion 
cyclotron waves and marginal sta- 
bility boundary for electrostatic 
ion cyclotron waves. 


the n = 1 wave as a function of the two remaining parameters, V^/a and T g /T^; 
the marginal stability contour for the k^ = 0 ion acoustic wave, [7], 

[9]; and the contour for y /w . = 1, the transition to the nonresonant 

*■ J 7 * ui^X pi 9 — - - ■* ■ — ■■ 

electron- ion beam acoustic instability. The shaded region denotes where the 
ion cyclotron wave is unstable but the ion acoustic wave is stable. The ion 
cyclotron wave is unstable to a smaller current over the range 0.05 < T /T . < 
8.5. Furthermore, it has a significant growth rate y /£!. = 0.1 (whill tfte 
ion acoustic wave is still stable) over the range 0.2? a ^ tVt^ < 2.5. Some 
modes of Tokamak operation may fall in this range. 
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3. 2 Electromagnetic Ion Cyclotron Waves. 


The full electromagnetic wave dispersion relation for spatially hono- 
geneous, drifting electron and ion Maxwellian distributions in a uniform 
magnetic field has been solved. In addition to the parameters above, the 
dispersion relation now depends upon 8^ = SimT^/Bp. For 3. large enough, 
electromagnetic ion cyclotron waves have a parallel phase velocity sufficient- 
ly slow that they can be destabilized by small electron drifts [5]. Fig. 5 
shows, for M/m = 1836 and T /T . = 1, contours of constant maximum growth rate 
in a 8 . , V /a parameter space, and the marginal stability contour y /Q. = 0, 
for the "electrostatic" ion cyclotron wave. This mode is electro- max 
statically polarized for 8. small, and develops small electromagnetic compo- 
nents as 8. increases. Finite 8. corrections are stabilizing. However, the 
electromagnetic ion cyclotron wave has a lower threshold for 8. > 10~ 3 . For 
Y /Q. = 10 -2 , real to - 0.75 ft . , and k^R. - 1.0 - 1.5, k |( R. < 1 0.1. As V 
decreases to where Y /Q. ^ 10 -3 , k»&. for maximum growth drops rapidly 

, . . „ , , r> max 1 X 1 

below 1.0 and k..R. increases. 

II i 

3 . 3 Nonlinear Considerations 

One- dimensional quasilinear electron plateau < f ormation can saturate the 
electrostatic ion cyclotron instability at e(j)/T 'v 10~ 2 , cf> = wave potential 
integrated over the unstable spectrum [4], [10]. Dum and Dupree [11], how- 
ever, have shown that if electron plateau formation is inhibited, ecp/T 'i 1/2, 
which agrees with Q machine results, [11], [12], [13], [14] where end 
effects may prevent plateau formation. We have observed [15] that resonance 
broadening due to large e<p/T permits ion cyclotron waves to interact with a 
significant portion of the electron distribution, leading to an effective 
electron- ion frequency v f - v ^ c ■ CT eJ; ) where A is the atomic number. 

In this limit, significant ion cyclotron anomalous resistance is possible. 

Experiments in the Princeton Model ST may operate in regimes where the 
electrostatic ion cyclotron wave is unstable, but where the ion acoustic wave 
may not be unstable. [16], [17], [18] It has been suggested that larger 8^ 
(higher density) may eliminate electrostatic current instabilities by reduc- 
ing the electron drift needed to produce the confining magnetic field. 

However, we must still contend with the electromagnetic ion cyclotron 
instability. Bers, et al. [19] have suggested that ion cycltron waves may 
account for the perpendicular ion heating observed in the TM-3 Tokamak. 

4. Turbulent Ion Heating 

Ion acoustic turbulent heating of a cold ion beam injected into a plasma 
with initial T /T. > 10 is observed in the UCLA DP (Double Plasma) device. 

A planar ion beam with drift energies exceeding T is accelerated from the 
driver plasma into the target plasma through a negatively biased grid. Fig. 

6 shows the background and beam ion energy distribution functions [23] for 
various axial distances from the grid, as measured by an electrostatic energy 
analyzer [23]. For the T = 1 eV case, the beam attenuation is due to charge 
exchange. The apparent change of background density is largely instrumental 
and arises from the increase in the volume from which background ions may 
arrive at the analyzer as it is pulled away from the injection point. For 
T =1 eV, and T. =0.2 eV, T /T . is below marginal stability for the ion-ion 
beam instability [24] and the ion beam propagates with only a slight (if any) 
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Fig. 6. Ion energy 
distribution vs. dis- 
tance. 


Fig. 7. Turbulent ion 
heating for various 
beam densities. 


Fig. 8. Turbulent ion 
heating for various 
beam energies. 


spreading in energy. For T^ = 3 eV, ion acoustic waves grow [2] , and the 
ion beam broadens and eventually forms a quasi-linear type of plateau. Beam 
ions heat both parallel and perpendicular to the beam direction; background 
ions also heat. 

In Fig. 7, the ion beam temperature vs. distance is displayed for 
several beam densities. The initial rate of ion temperature increase is 
approximately equal to the spatial linear growth rates of the ion acoustic 
waves. The dominant growing modes have wave numbers from 0.2 k to 0.6 k , 
and corresponding spatial growth rates of 0.01 k to 0.03 k . D The ion ° 
acoustic waves saturate at a rms density fluctuation level of 6n/n 0.01 to 
0.03. Ion heating continues after wave saturation, with the predominant 
beam kinetic to thermal energy transfer occurring in a distance 50A to 100A 
beyond the distance of maximum wave amplitude. In Fig. 8 the spatial growth 0 
of the beam temperature is displayed for various beam energies. The heating 
rates are reduced with increasing V /c as expected from the decrease in 

the ion-ion instability growth rate S with S V /c . 

beam s 

Directional probes show that when V > 2 c , waves grow obliquely to 
the beam, and that the wave energy densitils parallel and perpendicular to 
the beam are comparable for dense beams. Second harmonics of the maximally 
unstable waves are also observed; since k - 0.5 k_, these are probably gen- 
erated by resonant mode couplings between oblique waves. 

The initial exponential increase in ion beam temperature can be computed 

from quasi-linear theory. Taking the (V - _V ) 2 moment of the quasi-linear 

diffusion equation and assuming one dimensional, time independent conditions, 

we find that the beam temperature T^ satisfies 

b 

3T b /3x = mV I Vk ( ^ * 4 " w) (4) 

b b k 

where is the constant beam current, N is the wave action, and y is the 

growth rate. satisfies the approximate relation V 9N,/9 x - 2y N . 
Substituting this into (4) and integrating, we find approximately k k 
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1 


(5) 


V x) s n 


b b k 


7 V N. (0) (k *V - oj) 
y g k b 


exp[ 2k^x] 


where k. = y, /V and only the x-dependence of N, was considered, 
l 'k g k 

To estimate the ion acoustic rms saturation amplitude 6n/nQ, we utilize 
the results of strong turbulence theory (a la Dupree) by equating the initial 
linear growth rate y to the effective turbulent collision frequency = 

(k 2 D) 1 / 3 - kc^ (6n/nQ ) l ' 2 where D is the "average" velocity space diffusion 
coefficient. For y <\j k^V^ , and the observed kVk 1/20 and V /c g ^ 2 , we 
find 6n/n'Q 'v 10 2 . After turbulent wave saturation, the damping distance 
of the beam is Ax i V^/v ff = ^t/kc g (fin/n ^)” 1 2 . For the observed k n.0.2 k 
to 0.6 k^, we then have §x 30Xp t0 10 2 A^, ln reasona ble agreement with the 
experimental results; we conclude that the ion beam dissipation and thermali- 
zation are consistent with turbulent diffusion theories. 


5. Collisionless Electrostatic Shocks 

We have previously reported [21] data on steepening of ramp density 
waves, generated in DP devices, into collisionless electrostatic shocks 
(n^/n^ - 1.2) where n^ and n^ are the densities behind and ahead of the shock. 
We report here extensions to turbulent and to larger amplitude shocks [22]. 

The turbulent transition regions are generated by the reflected ion beam. 
Either ambient ion fluctuations or injected ion waves can be amplified in the 
ion reflection region to give turbulence, an effect considerably enhanced by 
the addition of small amounts of a light ion species (Figs. 9 and 10). 
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Fig. 9. Growth of natural ion Fig. 10. Amplification of an 

acoustic turbulence in a helium/argon injected wave train as it traverses 

mixture (n TI /n, = 3%). the shock (n„ /n. = 3%). 

He A He A 

When this occurs, the ion beam is scattered, resulting in a broadened distri- 
bution function. 

The large amplitude shocks are generated by a photoionized slab (n^ - 10 9 ) 
superimposed on a background of variable density, n^. At a ratio R (= n^/n^) 
above 30, the slab expands freely, with no shock formation. Steep shocks 
begin to appear around R = 10, and propagate with density jumps n./n^ - 2. 

The jump decreases as R decreases, eventually overlapping our DP data when 
R = 1.5 and n^/ng ~ 1.2. Typical profiles are seen in Fig. 11. The trailing 
wave train, typical at Mach 1.1 to 1.2, is heavily damped at higher Mach 
numbers, e.g., 1.5. The trailing wave structure is also sensitive to T /T., 
essentially vanishing when T /T. drops below 6. It is also very sensitive to 
the addition of light ions. Fig. 12 shows the effect of adding 3% helium to 
an argon plasma. 
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Fig. 11. 


Shock structure 


as a function of R = n /n_. 

R 0 



Fig. 12. Damping of trailing 
wave structure by light ions 
(He in A). An/n = 0.1. 


Fig. 13. Critical Mach number, M , 

for transition from small to large 

shocks , with M = V , , / (T /M ) 1 ' 2 . 

shock e A 


For a helium/argon mixture, we have studied the shock structure using 
cold fluid theory for the argon ions and a Vlasov treatment of the helium 
ions. The shock leading edge is simply described by the Sagdeev potential 
V(4>) , [25] where d 2 cj)/dx 2 = -dV/dcj), <j> = electrostatic potential. Above a 
critical Mach number, M , which depends on helium concentration and the ratio 
T e/ T i aS shown in Fi S* *£3, there exists a large amplitude shock which reflects 
almost all the helium ions and exhibits a double humped structure. Below 
the critical Mach number, M < M , most of the helium ions pass over the 


potential maximum which is thereby reduced approximately by the ratio 

(Below the lowest line, no shock exists in the He/A mixture.) This 

effect is present only if the helium ions are relatively cold (T /T > 

e He 


“h^V 
20 ). 
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